Recovery of chaotic tunneling due to destruction of dynamical localization by external 
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Quantum tunneling in the presence of chaos is analyzed, focusing especially on the interplay 
between quantum tunneling and dynamical localization. We observed flooding of potentially existing 
tunneling amplitude by adding noise to the chaotic sea to attenuate the destructive interference 
generating dynamical localization. This phenomenon is related to the nature of complex orbits 
describing tunneling between torus and chaotic regions. The tunneling rate is found to obey a 
perturbative scaling with noise intensity when the noise intensity is sufficiently small and then 
saturate in a large noise intensity regime. A relation between the tunneling rate and the localization 
length of the chaotic states is also demonstrated. It is shown that due to the competition between 
dynamical tunneling and dynamical localization, the tunneling rate is not a monotonically increasing 
function of Planck's constant. The above results are obtained for a system with a sharp border 
between torus and chaotic regions. The validity of the results for a system with a smoothed border 
is also explained. 
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I. INTRODUCTION 

Since the dawn of quantum theory, the tunneling effect 
has been recognized as a representative of genuine quan- 
tum phenomena. This does not necessarily imply that 
the tunneling effect has nothing to do with underlying 
classical dynamics. What recently attracts much interest 
concerns the fact that quantum tunneling in multidimen- 
sional systems or more specifically in non-integrable sys- 
tems strongly reflects the phase-space structures of the 
corresponding classical system 0, d, 0, HI, @, 0, @, 0, ED, 

EI. 

What is specific in generic nonintegrable systems which 
are neither completely integrable nor fully chaotic is that 
phase space is divided into infinitely many invariant com- 
ponents and they are generally intermingled with each 
other. The transition between different invariant com- 
ponents is forbidden in classical dynamics, and invariant 
components form dynamical barriers in phase space. On 
the other hand, in quantum mechanics, the transition 
through dynamical barriers thus formed may be allowed 
as a result of quantum effects. In particular, based on 
the notion of quasimodes [12], the quantum transition 
between congruent tori which are symmetrically formed 
in phase space is called dynamical tunneling in the liter- 
ature pjjj]. 

Dynamical tunneling happens in mixed-types phase 
space where invariant tori and chaotic regions coexist. 
This fact immediately invokes a question about the pos- 
sible roles of chaos in the process of dynamical tunnel- 
ing. The notion of chaos- assisted tunneling has been 
proposed to capture the situation where, besides quasi- 
doublet states localized on symmetrically formed congru- 
ent tori, chaotic states take part in the tunneling process 
0]. Semiclassical analysis in the time domain also re- 
veals that a bunch of tunneling trajectories associated 



with chaos, or more precisely the Julia set in complex 
plane, is involved in dynamical tunneling in mixed phase 
space 0, EH • What was commonly emphasized there was 
an aspect that chaos enhances dynamical tunneling. 

On the contrary, we here focus on another aspect of 
the role of chaos. In particular, it may happen that dy- 
namical tunneling is suppressed even in the presence of 
chaos. As explained below, there exist two competing 
effects behind the dynamical tunneling process, both of 
which originate from the nature of complex classical tra- 
jectories describing the tunneling transition from torus 
to chaos. 

As shown in Refs. 0, E3, exponentially many complex 
trajectories, not a single instanton path as in integrable 
tunneling, contribute to the semiclassical time evolution 
operator which allows us to evaluate the tunneling am- 
plitude from torus to chaotic regions. Such complex tra- 
jectories start from an initial state placed in the torus 
region and reach a final state in the chaotic sea. An im- 
portant fact is that in the first stage in their itinerary 
these complex orbits have non-small imaginary compo- 
nents and move in the complex plane, whereas the orbits 
are soon attracted by the real plane along stable mani- 
folds and then move very close to the real plane in the 
later stage. The rate of the attraction is exponentially 
fast, reflecting the fact that the behavior of the orbits is 
controlled by the stable manifolds of unstable periodic 
orbits in the real phase space 0, E3 • These observations 
imply that complex trajectories describing the tunneling 
transport have amphibious character: they act as tunnel- 
ing orbits when they stay in the torus region and behave 
as if they are almost real orbits when reaching the chaotic 
sea [15j. A well recognized aspect of dynamical tunnel- 
ing, that is, the enhancement of tunneling, is explained 
by the fact that the number of tunneling orbits is expo- 
nentially many 0, E3 • 
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However, their contributions to the tunneling ampli- 
tude are merely "potentially existing". This is because 
the time evolution in the later stage may involve oppo- 
site effects on the transport property, especially when 
dynamical localization gov erns the quantum dynamics in 
the chaotic region [I6LTI7L [H[ . The amphibious nature 
of complex orbits predicts that the tunneling penetra- 
tion through integrable barriers and the dynamics in the 
chaotic sea are not independent of each other, rather they 
must be closely related. Therefore, we expect that if 
the spreading of wave function is suppressed due to dy- 
namical localization in the chaotic sea, which appears 
as a result of destructive interference effects, the tun- 
neling transport is simultaneously suppressed. In order 
words, potentially existing tunneling orbits are exponen- 
tially many, but they interfere with each other to form 
localization. This is our working hypothesis and the pur- 
pose of this paper is just to verify that there exists strong 
interplay between dynamical tunneling and dynamical lo- 
calization. 

In order to confirm our hypothesis, that dynamical 
tunneling is significantly reduced as a consequence of dy- 
namical localization, our strategy is to apply noise to 
destroy the destructive interference generating dynami- 
cal localization. The idea is essentially the same as the 
one which has been used to see that the recovery of clas- 
sical diffusion occurs when external noise is applied to 
dynamically localized states [TH, [2(j. In the present case, 
to avoid noise- activated hopping from torus to chaos, we 
selectively add noise only to the chaotic region with the 
torus region being untouched. 

We should mention that the interplay between dy- 
namical tunneling and dynamical localization is crucial 
to understand the emergence of amphibious eigenstates 
[2lL [22| . In [111, [22[ they suggested the competition be- 
tween the effective Heisenberg time and the tunneling 
time is important to determine whether eigenfunctions 
show amphibious nature. The former time scale concerns 
dynamical localization and the latter dynamical tunnel- 
ing. Therefore, the flooding of eigenfunction found in 
|2lL [22I ] can be interpreted as a consequence of the flood- 
ing of potentially existing tunneling orbits. 

This paper is organized as follows. In Sec. [Til we 
present an area-preserving map that is designed to have 
phase space clearly separated into torus and chaotic re- 
gions. To verify our hypothesis stated above, a noise 
term is added. The noise term is characterized by two 
parameters: the strength e and the distance I from the 
border between torus and chaotic sea. In Sec. Ell we first 
investigate the case without noise as a reference for our 
subsequent argument about the case with noise. We espe- 
cially show two typical behaviors of tunneling transition 
and the sensitive dependence of the tunneling amplitude 
on the variation in the system parameters. In Sec. HV( we 
provide numerical evidence demonstrating that the intro- 
duction of the noise qualitatively changes the nature of 
tunneling. We define a tunneling decay rate 7, which 
will be used to quantify the rate of tunneling penetra- 



tion, and then show that it is strongly enhanced if noise 
is imposed. In Sec. El we examine how the noise strength 
affects the tunneling decay with the smallest /, i.e., the 
noise being applied to the whole chaotic sea. For small 
e, 7 increases pert ur bat ively as e 2 increases. For large 
e, 7 saturates to a certain "classical value" and the pa- 
rameter sensitivity of 7 disappears. In Sec. IVI[ we focus 
on the /-dependence on 7, where we show that a char- 
acteristic scale of I coincides with the localization length 
of the chaotic sea. An interesting finding is that there 
exists a region of I in which the tunneling rate decreases 
even though Planck's constant becomes large. This is ex- 
actly because of the strong interplay between dynamical 
tunneling and dynamical localization. In Sec. rVTI[ we ex- 
amine whether or not our result is valid for systems with 
a smoothed border between torus and chaotic regions. 
We show that almost all results obtained in cases with 
clearly separated phase space hold even in such generic 
situations. Finally, in Sec. I Villi we discuss our results 
and summarize the paper. 



II. DESIGNING PHASE-SPACE FOR CHAOTIC 
TUNNELING 

Classical phase space in a mixed system is very com- 
plex in general; besides chaotic components, a variety of 
invariant components such as Kolmogorov-Arnold-Moser 
(KAM) circles, islands of stability, or cantori live together 
in a single phase space. Great difficulties arise not only in 
classical but also in quantum dynamics in mixed systems 
due to the immense complexity of phase space. 

In order to avoid complicated effects due to the mixture 
of different types of invariant structures, one strategy 
is to make the phase space "clean." Designing sharply 
divided phase space in a class of piecewise linear map is 
in this spirit [23[, and the recent finding of mushroom 
billiards also provides a good testing ground for a precise 
understanding of mixed systems [241 ] . 



A. Model system 

We here present our system whose classical phase space 
is clearly separated into torus and chaotic regions. Let 
us consider the one-dimensional kicked rotor described 
by the Hamiltonian 

00 

H(p,q,t)=T(p) + V(q) S(t-n), (1) 

n— — oo 

where p and q are the momentum and the position (an- 
gle) of rotor, respectively. The corresponding classical 
dynamics is reduced to the area preserving map 

F • ( Pn+l \ _ ( Pn~ V'(q n ) \ 

'Un+l J ~ Un + T'(p n+1 ) J' ^ 
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FIG. 1: (Color online) Phase space portraits of the classical 
map JTJ) with T(p) — To(p) (Eq. ©). The parameters are set 
to be d = 0.3, u = 0.6418 • • • , k = 2 and (a) s = 0.5, (b) 
s = 2, and (c) s = 4. The smoothing parameter (3 = 100 is 
chosen so that the border between torus and chaotic regions 
is clear. Note that the border is not sharp in the strict sense 
(see text). A magnification of the area around the border in 
the case of s — 4 is shown in (c'). For small s, there are big 
islands in the chaotic sea and also around the border. On the 
other hand, if s is sufficiently large, e.g., s = 4, these islands 
become so small that the influence of these almost invisible 
islands on dynamical tunneling is negligible. 

To realize the phase space which is clearly separated 
into regular and chaotic regions, we take the kinetic and 
potential terms in the following way: for the kinetic term, 
we have 

To(p) = \s(p ~ d) 2 p (p - d) + u(p - d), (3) 

where 0p(x) denotes a smoothed step function with a 
parameter /3(> 0): 

0p(x) = ±{l + tanh(l3x)}, (4) 

e.g., 0p{x) — > as x — > — oo and 0p(x) — > 1 as x — > +00. 
lim^^oo 6p(x) gives the Heaviside step function. For p <C 
d, To(p) ~ uo (p — d) is a kinetic term for a linear rotor 
with a rotation number uo. On the other hand, above the 
"border" p = d in the momentum space, the quadratic 
term s(p—d) 2 /2 appears in the kinetic term. Throughout 
this paper we employ the following potential term 

V{q) = ±coB(2*q). (5) 
We therefore have, for p <C d, 
k 

H « u(p- d) + cos(27r<2) S(t-n) (6) 

n— — oo 

which is completely integrable, and, for p^> d, 

1 k °° 

- s (p-d) 2 + u(p-d) + — cos(27rg) ^ 5(t-n), 

n— — 00 

(7) 

which is reduced to the standard mapping with a nonlin- 
ear ity parameter sk under an appropriate scaling of the 



Hamiltonian and a shift in the momentum space. Fig- 
ure [T] shows examples of the classical phase space with 
several different values of s.In the following, we will ex- 
amine the case of s = 4 where phase space is clearly 
separated into torus and chaotic regions. As shown in 
Fig. [1] (c), when s = 4 the islands around the border be- 
tween torus and the chaotic regions are almost invisible. 
However, it should be noted that the border is not as 
sharp as in a previously studied class of piecewise linear 
map [23] or mushroom billiards (24|. One can find tiny 
island chains in much more magnified figures. 

B. Introduction of the effect of noise 

As is well known, dynamical localization occurs as 
a result of destructive interference in the chaotic sea 
[HL [13, EH- Conversely, once quantum interference is 
destroyed in some way, classical diffusion is recovered 
pH [20| . In a similar way, we will verify that chaotic 
tunneling is suppressed by dynamical localization in the 
chaotic sea by examining whether or not potentially ex- 
isting chaotic tunneling is recovered if destructive inter- 
ference creating dynamical localization is attenuated. 

A simple way to destroy quantum interference is, as 
done in numerical tests to see the recovery of classical 
diffusion [l9|, [2O], to apply noise to chaotic seas. How- 
ever, the introduction of noise needs to be done carefully 
and should only be limited to chaotic regions, otherwise 
the tunneling process will be contaminated by the ther- 
mal hopping between the torus and the chaotic seas. In 
order to avoid the transition from the torus region into 
chaotic sea, noise will be applied only to the region that 
is far from the torus, say p > L (> d). We employ the 
following modification to the kinetic term to realize such 
momentum-dependent noise: 

T e (p, n) = T (p) + eUP ~ L)0p(p - L), (8) 

where e is the noise strength and £ n is a stochastic vari- 
able. We assume that £ n obeys a Gaussian distribution: 

(£n> = 0, (Zntim) = ^mn, (9) 

where (•) denotes the ensemble average. 

C. Quantum map 

The quantum dynamics associated with the classical 
map Eq. ([2]) is presented following the standard quan- 
tization procedure of quantum maps. Due to the pe- 
riodicity of the map Eq. ([2]) in the periodic bound- 
ary condition (with a period 1) is imposed on the wave 
function in the position representation. We also impose 
the periodic boundary condition in the momentum space 
— \ <P< \- As a result of these boundary conditions, 
the effective Planck's constant is thus quantized as 

W 
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FIG. 2: Probability distribution of initial wave function in 
the momentum space (upper) and its schematic illustration in 
phase space (lower). We impose noise on the shaded domain. 

where N is the dimension of the Hilbert space of the 
system. Furthermore, in order to avoid the finite size 
effect in the diffusion process, an absorbing boundary is 
set at p = p C utoff- We deliberately choose p cu toff large 
enough to ensure the convergence of the characteristics 
of dynamical tunneling from torus to chaotic sea [25[ . A 
unit time evolution of the state vector \ip n ) is described 
by the quantum map 

|V„+i) = u\ip n ), (ii) 

where U is the time evolution operator for a unit time 
interval 

U = exp (-^T(p)) exp ("^(<z)) • (12) 

III. TUNNELING DOUBLET, SUPPRESSION 
OF DECAY, AND DYNAMICAL LOCALIZATION 
IN THE SYSTEM WITHOUT NOISE 

Before studying the case where the noise term is 
present, we show that the tunneling process from the 
torus to the chaotic region is sensitive to the change in 
system parameters, even though the classical phase space 
of our system has carefully been designed to remove res- 
onance islands. The result tells us that the quantum in- 
terference effect is very relevant in the issue of dynamical 
tunneling, and it strongly affects the transition process 
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FIG. 3: (Color online) An example of the off-resonant case 
(s 4.2566 •••) with h 1/20, W = 40 and p cuto ff 18. 
Other parameters are the same as Fig. [T] At n — almost all 
of the probability is localized in the torus region, (a) P* orus as 
a function of time step n. During the time evolution, P* orus 
keeps almost unity value. In the inset, we magnify the tiny 
oscillation of P^ orus . (b) Snapshots of the momentum distri- 
bution \ipn(p)\ 2 - The border between the torus and chaotic 
regions is at d = 0.3. The diffusion in the chaotic sea is sup- 
pressed by dynamical localization, so the probability distribu- 
tion tends to a stationary distribution with an exponentially 
decaying tail, which is characteristic of dynamical localization 
in the chaotic sea. 

Hereafter, the time evolution of tunneling is quantita- 
tively monitored using 

rPb 

Pn ms = / \MP)\ dp, (13) 

where ij) n {p) is the wave function in the momentum rep- 
resentation at the n- th step, and parameters p a and p^ 
are appropriately chosen to enclose the torus region. We 
use p a = — W/2 and p = d in the following [2r| . 

Throughout this paper, the initial wave function is 
taken as an eigenstate of the linear rotor Eq. ([6]). As 
explained in Appendix [A[ the initial wave function 
ip torus (q; p c ) in the position representation is specified by 
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FIG. 4: (Color online) An example of the resonant case. Ex- 
cept for s( = 4.2178 • • • ), the values of the parameters are the 
same as in Fig. [3] (a) P^ orus as a function of time step n. A 
large oscillation ranging from 0.1 to 1.0 is observed, (b) Snap- 
shots of the momentum distribution \^n(p)\ 2 at n — 200, 1800 
and 3500, each of which corresponds to an extremum in the 
curve shown in (a). Dynamical localization occurs in the 
chaotic region, as it does for the off-resonant case. In the 
inset, tunneling oscillation between the torus region p < d 
and a part of the chaotic sea d < p < 5 are magnified. 



Pc where p c represents the "average" of the momentum 
of the torus [ItJ. The initial state is placed just below the 
border between the torus and the chaotic sea, i.e. p c = 
(see Fig.[2j). This is because the tunneling amplitude be- 
comes too small compared to the numerical precision if 
the initial state is placed too far from the border. 

As mentioned above, even with a small change in pa- 
rameters, P^ orus changes rather drastically and there ap- 
pears a variety of oscillatory patterns in the plot of P* orus . 
In the following, we show two typical examples. As seen 
below, in spite of only slightly different values of s, oscil- 
latory patterns are entirely different. Note that variation 
in s does not have any effect on the integrable region. It 
only affects chaotic regions. The classical phase space 
portrait is therefore almost unchanged under a small 
change in s. 




FIG. 5: The oscillation amplitude of P n orus as a function of 
s. 
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FIG. 6: Quasienergies of torus and chaotic states (upper) and 
the magnification of oscillation amplitude of P^ orus in Fig.[5]as 
a function s (lower). Only the quasienergies of the eigenstates 
whose momentum expectation values fall within the range 
—0.2 < p < 6.0 are shown [29], because the chaotic states 
which have large contribution to the oscillation are expected 
to be near the torus state. A line running in the horizontal 
direction represents an eigenenergy associated with a torus 
eigenstate. It forms many avoided crossings with energy levels 
corresponding to chaotic states. Note that peaks of oscillation 
amplitude appear when avoided crossings occur. 



First we depict the time evolution of P* orus in the case 
of s = 4.2566 • • • (see Fig. [3]). We hereafter call this the 
off-resonant case [28j. P* orus keeps almost unity value as 
shown in Fig. 02(a), which means that most of the proba- 
bility amplitude is confined within the torus region. On 
the other hand, as shown in Fig. [3] (b), the snapshots of 
the log-scaled probability distribution in the momentum 
representation, once the tunneling transition from the 
torus to chaotic region occurs, the wave function starts 
to diffuse into the chaotic sea. However, the diffusion 
is gradually saturated and the wave function asymptot- 
ically approaches a stationary shape. This is a result of 
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dynamical localization. In other words, almost all of the 
probability amplitude that tunnels out of the torus region 
is confined in a region d < p < d + £. Here £ represents 
the localization length in the chaotic sea 

Next, we slightly change the parameter to s = 
4.2178- •• (see Fig. Hj). In contrast with the previous 
case shown in Fig. EKa), P* orus exhibits a periodic oscil- 
lation (see Fig. Ufa)). In accordance with the oscillation 
of P* orus , as shown in Fig.|H(b), the momentum distribu- 
tion also exhibits a recurrent pattern. Initially, the wave 
function is localized within the torus region. As time 
elapses, most of the probability amplitude escapes to the 
chaotic region. After that, the wave function comes back 
to the torus region. Such a recurrent oscillation strongly 
suggests the emergence of a tunneling doublet composed 
by a torus and a chaotic state. We call this type of be- 
havior the resonant case [281 ] . 

To see large fluctuation with the variation in s more 
explicitly, we present in Fig. [5] the amplitude of the time- 
dependent oscillation of P* orus as a function of s. Such 
an erratic oscillation pattern reminds us of chaos-assisted 
tunneling, where the tunnel splitting between symmetry- 
related torus states exhibits strong fluctuation with vari- 
ation in parameters [l|, [3[ . In the case of chaos-assisted 
tunneling, the fluctuation is induced by avoided crossings 
that involve the torus and the chaotic states. 

Now we explain the origin of the sensitivity of the tun- 
neling to the variation in s in the present case. This is 
also due to the sensitivity of accidental degeneracies be- 
tween torus and chaotic states. As evidence, we compare 
the quasienergies and the oscillation amplitude of P* orus 
in Fig. [6l Since s mainly affects the dynamics in the 
chaotic sea, the quasienergies of the torus states keep al- 
most constant and form avoided crossings with many en- 
ergy levels corresponding to chaotic states. Around these 
avoided crossings, coherent oscillation between torus and 
chaotic states occurs. In Fig. [6j we clearly see the cor- 
respondence between peaks of the oscillation amplitude 
and positions of avoided crossings. Note that the period 
of the tunneling oscillation is well described by the for- 
mula T osc = 27rft/A, where A is the gap of quasienergies. 

We remark on the influence of the absorbing bound- 
ary. The distance between the positions of the border d 
and the absorbing boundary p cu toff is much larger than 
the localization length, e.g., p C utoff — d ^> £. The prob- 
ability amplitude which reaches the absorbing boundary 
is exponentially small, even if almost all of the proba- 
bility amplitude temporarily gets out of the torus region 
in the resonant case (see Fig. IU[b)). Hence the influence 
from the absorbing boundary is virtually negligible. In 
this sense, the wave function is approximately bounded 
around the torus region or the nearby chaotic sea. 

We summarize our findings in this section. In the ab- 
sence of noise, the behavior of the wave function launched 
from the torus region sensitively depends on the parame- 
ter s while classical dynamics does not. This is explained 
by the presence of the near-degeneracies of quasienergies. 
In particular, we discovered a new kind of tunneling dou- 
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FIG. 7: (Color online) The off-resonant case under the pres- 
ence of noise with L — d + 0.2. The other parameters are 
the same as in Fig. [3] (a) P* orus exhibits exponential decay, 
where the exponent increases as e increases, (b) We depict 
snapshots of the probability distribution at e = 0.005. For 
large n(> 100), the effect of noise becomes significant (cf. 

Fig. Eft)). 



blets that couples a torus state with a chaotic state that 
exhibits dynamical localization. 



IV. EXPONENTIAL DECAY DUE TO THE 
BREAKDOWN OF DYNAMICAL 
LOCALIZATION 

As shown in the previous section, the tail of the wave 
function, which is initially localized and exponentially 
decaying towards the chaotic sea, gradually spreads over 
the chaotic sea as time proceeds and dynamical localiza- 
tion controls the nature of the wave function there. In 
other words, dynamical tunneling necessarily accompa- 
nies dynamical localization in the chaotic sea. 

To examine the influence of dynamical localization on 
the tunneling process further, we will add external noise 
to the chaotic region to attenuate dynamical localization. 
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FIG. 8: (Color online) The resonant case with noise (L = 
d + 0.2). The other parameters are the same as in Fig. [7] 
(a) P* orus exhibits oscillatory decay for e = 5 x 10" 4 . To 
introduce a decay rate, we fit an exponential curve to the 
oscillatory decay for short time steps (e.g. n < 10 4 ). A 
dashed line in the inset depicts the result of the fitting to 
the oscillatory decay. As the effect of noise increases, the 
resonant oscillation disappears and P^ orus shows exponential 
decay for longer time steps, where the decay rate slows down 
in contrast with the off-resonant case. The decay rate for the 
strongest noise (e.g. e = 0.005) is almost the same as that for 
the off-resonant case, (see Fig. [7(a)). (b) Snapshots of the 
time evolution at e = 5 x 10 -4 are presented. Each time step 
corresponds to the time step in Fig. HJb). 



As is explained in Section IIIBj the kinetic term of our 
kicked Hamiltonian is replaced with T e (p, n) (Eq. (|HJ)), 
where noise is applied only to the region p > L. To 
avoid contaminating the tunneling process with "thermal 
agitation" , we impose the condition L > d. Furthermore, 
in this section, we make L as small as possible to observe 
the case where dynamical localization is attenuated in 
the whole chaotic sea. In the following, we examine the 
effect of noise for both resonant and off-resonant cases. 

First, for the off-resonant case (s = 4.2566 •••), we 
can explicitly see from Fig. 0(a) that unlimited diffusion 
does take place [30|, which is in sharp contrast to Fig.[3j 



More importantly P* orus decays exponentially, and its de- 
cay rate depends on the intensity of the noise. As seen in 
Fig. [71(a), the decay rate gets larger with the increase in 
the noise intensity. The origin of the exponential decay of 
ptorus can 1^ attributed to the tunneling effect, since we 
apply noise so that it does not affect the classical dynam- 
ics in the torus region. In other words, the exponential 
decay has nothing to do with thermal hopping from the 
torus into the chaotic regions. Typical snapshots of the 
probability distribution are shown in Fig. [3(b), where the 
effect of noise becomes significant at larger n. 

Thus we show for the first time that chaotic tunnel- 
ing from a torus into chaotic regions is enhanced drasti- 
cally even though the noise is added only to the chaotic 
sea. This means that the recovery of diffusion in the 
chaotic region or the disappearance of dynamical local- 
ization gives rise to drastic enhancement of dynamical 
tunneling, or alternatively stated, dynamical localization 
suppresses dynamical tunneling. 

Second, we examine the resonant case (s = 4.2178 • • • ), 
where P* orus periodically oscillates in time. In the pres- 
ence of noise, the behavior of P* orus is sensitive to e (see 
Fig. [81(a)). P^ orus decays much faster than in the off- 
resonant case (see Fig. [71(a)). This is because dynamical 
tunneling, whose restoration is due to the presence of 
noise, as in the off-resonant case, is pumped by the reso- 
nance, which remains because the noise intensity is still 
small enough. Snapshots of probability distribution are 
shown in Fig. [51(b). We also show a magnified view of 
Husimi functions in Fig. [9] to explain the pumping. For 
larger e, P* orus exhibits exponential decay. Furthermore, 
as e increases, the decay of P* orus becomes slower, in con- 
trast with the off-resonant case. This suggests that the 
effect of resonance becomes weaker as the noise intensity 
becomes larger. 

These observations motivate us to introduce a de- 
cay rate 7 for P* orus to characterize the restoration of 
dynamical tunneling, at least in an early stage, e.g., 
n < 10 4 - 10 6 steps in Figs. and We should re- 
mark, however, that the fitting is rather subtle for weak 
e in the resonant case, where the "early stage" becomes 
rather short, e.g., n < 10 4 in the example shown in Fig. [SI 



V. NOISE STRENGTH DEPENDENCE OF 
CHAOTIC TUNNELING 

In this section, we examine the dependence of the tun- 
neling decay rate 7, introduced in the previous section, 
on the noise strength e. We first recall the studies on 
the effect of noise on dynamical localization [l9|, [20[ , in 
which it was reported that the diffusion rate of quan- 
tum systems when noise is imposed depends on the noise 
strength. More precisely, if noise is weak enough, the 
behavior of the diffusion rate as a function of the noise 
strength is understood within a perturbation theory: the 
diffusion rate scales with the square of noise strength. In 
contrast, if noise is stronger than a certain threshold, the 
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FIG. 9: (Color online) Time evolution of Husimi functions 
(in a logarithmic scale) in the resonant case for e — (left 
column) and e = 0.0005 (right column) at n = 200, 1800, 3500 
and 5000 (from top to bottom). The corresponding phase 
space portrait is superposed. The other parameters are the 
same as in Figs. [4] and [8] For e = 0, there occurs coherent 
tunneling between the torus region and the chaotic sea. At 
n = 1800 and n — 5000, the system completely tunnels out to 
a localized region in the chaotic sea. The system periodically 
comes back to the torus region (e.g., n — 3500). Hence the 
"pumping out" and "pumping back" of wavepacket between 
the torus region and chaotic seas cancels out in the absence of 
noise in the long time average. The pumping still remains for 
e — 0.0005, though the cancellation is broken. On a short time 
scale, the effect of noise is negligible (n = 200). At n = 1800, 
the noise hinders the tunneling from the torus partially, and 
at the same time, promotes the diffusion at the chaotic sea. 
The discrepancy due to the noise is more evident at n = 3500. 
Although Pn=3500 is far larger than Pn=iIoo (see the inset of 
Fig- 0(a)), the recurrence is not quite perfect. In the chaotic 
sea, there remains a diffusive component, which is produced 
by the tunneling in the previous steps and fails to back to the 
torus region. This explains why the pumping that occurs in 
the resonant case promotes the tunneling decay (see, Figs. [7] 
(a) and 0(a)). 



diffusion rate saturates to a classical value. We here ex- 
pect that a similar dependence of tunneling on the noise 
intensity would be observed. Figure [10] actually verifies 
this prediction. Below we will explain how 7 depends on 
e in more detail. 

When noise is weak (e/h <C 1), the decay rate 7 is 
proportional to (e/h) 2 . This behavior can be, although 
not explicitly shown here, reproduced by performing a 
perturbation expansion with respect to a small parame- 
ter e. Qualitatively, exponentially decaying behavior is 
observed as a result of destruction of dynamical localiza- 
tion, therefore, as the degree of destruction is increased, 
the decay rate 7 grows. Furthermore, 7 strongly depends 
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FIG. 10: (Color online) (a) The decay rate 7 as a func- 
tion of e/h for different values of the parameter s = si(= 
4.175 • • • ),s 2 (= 4.2178 • • • ) and s 3 (= 4.2566 • • • ). The for- 
mer two correspond to resonant cases and the latter to an 
off-resonant case. The other parameters are the same as in 
Fig. When the effect of noise is weak, the tunneling decay 
rate obeys a perturbative behavior 7 ~ e 2 (dotted line). 7 is 
saturated in the strong intensity regime (e/h > 1). (b) The 
decay rate 7 as a function of s for different e/h. Arrows indi- 
cate Sj's. The locations of the peaks for small e coincide with 
the resonances of the unperturbed system (see also Fig. [5}. 



on the parameter s, and in particular, the locations of the 
peaks for finite e coincide with the resonances of the un- 
perturbed system (e = 0), which is shown in Fig. [10] (b). 
This means that the effect of resonance persists in the 
perturbative regime. 

When e lies in an intermediate region, the decay rate 
7 is still influenced by the effect of noise. In the off- 
resonant case, 7 obeys the perturbative behavior over a 
rather wide range of e. On the other hand, in the res- 
onant case, the perturbative argument is not applicable 
in the intermediate region, and the rate of growth of 7 
with e depends on s. Two typical examples are shown 
in Fig. [10] (a) ( 10~ 2 < e/h < 1 ). With s = si, the 
rate of growth of 7 with e 2 temporarily decreases and 7 
converges to the value for the off-resonant case (s = S3) 
as e is increased. In another case, s = S2 in Fig. [10] (a), 7 
temporarily decreases even with increase in e. A possible 
explanation is that the effect of resonance is attenuated 
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FIG. 11: (Color online) (a) The logarithm of the tunneling 
rate 7 as a function of / = L — d for several values of s. 
(b) The average and variance of log 10 7. A theoretical line 
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(see text) is drawn as a dashed line. The 



parameters are h— 1/20 and e/h = 5. 



by moderate noise. 

We summarize the results for weak and intermediate 
noise intensity regions: an overall trend is that 7 in- 
creases on average, whereas the fluctuation of 7 decreases 
with increase in e. This suggests that, in the weak and 
intermediate regions, there exists a competition of two 
different interference effects; one is a constructive inter- 
ference induced by resonances, and the other a destruc- 
tive interference associated with dynamical localization. 
The former and the latter correspond to the changes of 
the fluctuation and the average of 7, respectively. 

For large e, it is natural to expect that these interfer- 
ence effects become totally incoherent. Indeed, we nu- 
merically find that 7 is almost independent of s, which 
is shown in Fig [lOl What is most important is that 7 
seems to converge to a certain value. The existence of 
the plateau regime (in the noise strength dependence of 
the tunneling rate) strongly suggests a "classical value" 
of the tunneling rate 7 in the incoherent limit. If this is 
the case, it must be evaluated based on complex classical 
orbits. 



VI. PARTIAL DESTRUCTION OF 
DYNAMICAL LOCALIZATION 

In the cases considered so far, noise was applied to the 
whole chaotic sea, that is L « d. In this section, we 
see how the decay rate 7 behaves and how the localiza- 
tion length £ depends on 7 when noise is applied only 
to the limited chaotic region which is distant from the 
torus region, that is L > d. In the following we denote 
the width of the unperturbed chaotic region, that is the 
region without noise, by I = L — d as shown in Fig. [2] 
schematically. 

For I > 0, since dynamical localization remains in the 
unperturbed chaotic region d < p < L, we can expect 
that dynamical tunneling is suppressed partially by the 
remaining localization. On the other hand, in the per- 
turbed region p > L, we can expect that noise destroys 
dynamical localization to restore the diffusion process for 
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FIG. 12: (Color online) The scaled exponent of the decay rate 
(log(7/7o)) s as a function of the scaled width Z/£. The dotted 
straight line is the relation (|14[) with a = 2. 
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FIG. 13: (Color online) Ensemble averaged decay rates 
(log 10 7)s for h = 1/5 (square) and 1/20 (circle). The the- 
oretical lines 7 ~ e~ 21 ^ are also drawn. The undershoot in 
the plot of h = 1/5 is an artifact of numerical fitting to the 
data from the limited time evolution, e.g. n < 10 7 . When 
the width I is large (I > 1), the decay rates with a smaller 
Planck's constant (h — 1/20) become larger than those with 
a larger Planck's constant. This is because the localization 
length £ increases as the Planck's constant decreases. 



the same reason as discussed above. The latter is again 
supposed to help the restoration of the chaotic tunneling. 
In order to avoid complications, we here restrict ourselves 
to the case where the noise intensity is sufficiently large, 
which corresponds to the case with large e examined in 
the previous section. 

As is shown in Fig. [TTJ we found that 7 tends to de- 
crease exponentially as I increases. Furthermore, when I 
is larger than a certain value, 7 strongly fluctuates with 
respect to s. An interpretation for the presence of large 
fluctuation is that the resonance previously observed in 
the perturbative regime (i.e., smaller e) with I « is re- 
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covered. It should also be noted that 7 is not monotonic 
with I. We suppose that this is also due to the recovery 
of the resonance, whose quasienergy gap is modulated by 
I. 

In the following, we discuss a gross feature of the 7- 
dependence on I through an ensemble average with re- 
spect to s within a narrow interval As = 0.05. Since 
the fluctuation of 7 is rather large, we take the average 
(In 7) s , instead of {j) s - As shown in Fig. [TTTb). the av- 
erage (In 7)3 decreases monotonically as a function of I. 
This is consistent with exponential dependence of 7 on 
I. On the other hand, the fluctuation of In 7 gradually 
becomes larger as I increases. This is due to the presence 
of resonances as explained above. 

The exponentially decaying behavior of 7 as a function 
of I suggests the presence of a characteristic momentum 
scale. It is simplest to hypothesize the following depen- 
dence of 7 on the localization length £ of the unperturbed 
chaotic region. 

(In 7) a ~ In 70 - (14) 

where 70 and a are constants. 

To examine the validity of this hypothesis, we first de- 
termine the localization length £ by fitting the slope of 
the probability distribution numerically. The value of a 
can be roughly estimated using the following perturba- 
tive argument. Let a wavepacket evolve from the torus 
without noise. If we wait long enough, the wavepacket 
becomes a superposition of the torus state and a chaotic 
state. The latter exhibits dynamical localization in the 
chaotic sea (see, Sec. IIII|) . In this case, the envelope of 
the wave function in the momentum representation takes 
the exponential form ~ e~ v '^ in the chaotic sea. Once 
we impose noise at p > L, tunneling penetration starts. 
We assume that the wave function for p < L is almost 
unchanged in spite of the presence of noise. Then the 
rate of the tunneling would be proportional to \e~ l ^\ 2 . 
This implies a = 2. A test of our argument for h = 1/20 
is shown in Fig. HITb). The result suggests a should be 
slightly larger than 2. 

In Fig. [12j the validity of Eq. ([T4|) for various values of 
h is examined. If a in Eq. ([T4]) is assumed to be indepen- 
dent of /i, 7 depends on h only through 70 and £. Hence 
we check how the scaled decay rate (111(7/70) ) s depends 
on the scaled width of the unperturbed chaotic region 
Z/£. The results in Fig. [12] roughly validate our assumed 
scaling relation Eq. ([T4]) with a slightly less than 2. We 
thus conclude that the localization length is a crucial fac- 
tor for determining 7. In other words, chaotic tunneling 
is shown to be strongly correlated with dynamical local- 
ization in such a way. 

We discuss an implication of the scaling relation ([T4|) 
for ^-dependence on 7. On one hand, 70 decreases as 
h —> 0. On the other hand, the localization length 
£ increases as h — > [l7| • Hence, h dependence of 
7 7 exp(— al/£) for a given value of I is not mono- 
tonic, provided I > 0. An example is shown in Fig. [131 
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FIG. 14: (Color online) Phase space portraits for small ft 
(= 5), where the border between torus and chaotic regions is 
blurred, and island chains appear around the border. Other 
parameters are given as d — 0.45 and s — 4. A full view and 
its magnification are shown in (a) and (b), respectively. In (a), 
we depict a Lagrangian manifold corresponding to the initial 
wave function ip torus (q; p c ) with p c = by a thick curve. Also, 
a thick horizontal line at p — d indicates the discontinuity of 
the smoothed step function that appears in the kinetic term 
To(p). When ft is large, the line p — d touches the border (see 
Fig. [TJ. On the other hand, for small ft, the line p — d and 
the boundary are detached [3^]. (b) Magnifications around 
the border without noise (e = 0.0, L = d + 0.2) and with 
noise (e = 0.005, L — d-\- 0.2). Note that there are almost no 
effects of noise around the border. 



where I vs. (111(7) ) a is plotted for two values of h. For 
small /, (111(7) ) a becomes smaller as h becomes smaller. 
This is explained by ^-dependence of 70. For large Z, 
the opposite occurs. Namely, the tunneling rate grows 
even though the Planck's constant h is reduced. Hence 
we discovered that the Planck's constant dependence of 
the tunneling rate is not monotonic when there partially 
remains dynamical localization in the chaotic sea. This 
is quite a unique feature and another piece of evidence 
for the strong correlation between dynamical localization 
and chaotic tunneling. 



VII. THE CASE WITH SMOOTHED BORDER 

Up to now, we have studied the system whose phase 
space is sharply separated between torus and chaos, and 
the degree of sharpness has been controlled by the pa- 
rameter ft. This is motivated, as stated in section [Til 
to remove complicated effects caused by the presence of 
complicated invariant structures in generic phase space. 
However, as a price we have to pay, the sharpness of 
the border between torus and chaotic regions may induce 
anomalous effects in quantum mechanics. 
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FIG. 15: (Color online) Dependence of amplitude of tunneling 
oscillation and tunneling decay rate on parameter s in the case 
of blurred border (a) amplitude of tunneling oscillation, (b) 
tunneling decay rate 7. The parameters of the underlying 
classical systems are the same as in Fig 1141 Here we take h = 
1/20. In (b), the strengths of e are the same as in Fig. [10] (b). 



What is important to realize a sharp border is to take 
the kinetic term Tq(p) with discontinuity (Eq. (pJJ)). Be- 
cause of the nature of the 0p(x) presented in (|4]), the 
second derivative T^(p) is nearly discontinuous. Espe- 
cially, in the limit (3 — > 00 Tq(p) is strictly sharp and thus 
the free evolution gives rise to the so-called Gibbs phe- 
nomenon in the Fourier transformation of the wave func- 
tion in the momentum representation. Namely, the wave 
function is strongly diffracted (32[. Furthermore, since 
the discontinuity of T$(p) lies near the border through 
which dynamical tunneling proceeds, diffraction can have 
strong influence on the tunneling process. In our calcula- 
tions, we have used a large, yet finite (3 (say, 100 in Fig. [I]) 
, but the observation of Husimi functions in the first sev- 
eral steps clearly reveals that there remains a signature 
of diffraction even in case (3 = 100 (33[. Therefore we 
here examine whether or not what has been claimed in 
the previous sections is valid even in small (3 cases, where 
the signature of diffraction is expected to be negligible. 

In the case of large (3 shown in previous sections, a 
torus state and a localized state in the chaotic sea play 
a central role. As (3 becomes small, resonance islands 
appear around the border and the island states associated 
with them need to be taken into account in addition to 
torus and chaotic states 0, 0, • This may invoke two 
kinds of "tunneling channels" (1) a torus to chaotic sea; 
(2) a torus to an island, first, and subsequently the island 
to chaotic sea. 

Our scenario for the numerical study of small (3 cases 
is almost the same as before: launch a wave function 
from an approximate torus state ip torus (q; p c ) and moni- 
tor the time evolution of P* orus . We add noise in a similar 
manner as before to destroy the quantum coherence pro- 
ducing dynamical localization. 

For (3 = 5, as shown in Fig. [HJ there emerges nonlin- 



ear resonances in the torus region and small islands in the 
chaotic sea as expected. It is important to note that small 
islands thus created are almost unchanged even though 
noise is applied in the region p > d. An example of a 
phase portrait is presented in Fig. [141(b). This suggests 
that the changes of the small islands at different noise 
strengths do not have any influence on the tunneling pro- 
cess. 

In the absence of noise, whether the quasienergy reso- 
nance occurs or not, the wave function is almost bounded 
in the torus region and only a small fraction of tunneling 
tail penetrates the dynamical barrier. Thus, the behavior 
of P* orus is similar as that shown in Figs. [3] and [4] As is 
seen in Fig. [15] (a), the amplitude of the tunneling oscil- 
lation is also sensitive to the variation in the parameter 
s. 

Note that the sensitivity with respect to s is not re- 
duced even with smaller (3. The origin of the sensi- 
tivity should be attributed to the resonance among ap- 
proximate stationary states, as explained in Section [V] 
We comment on the contribution from the island states, 
which was not taken into account in the explanation in 
Section [V] This is expected to be small, since the varia- 
tion in s influences only on the region far apart from the 
torus region, and the quasienergies of island states do 
not depend on s. If a quasienergy resonance is sensitive 
with respect to s, it must involve a localized state in the 
chaotic sea. Hence we conclude such a resonance is the 
origin of the strong fluctuation. 

As shown in Fig. [15] (b), we have confirmed that chaotic 
tunneling is restored even for small /?, and that the be- 
havior of the decay rate 7 also persists. We note that 7 
for small (3 needs to be deliberately determined from the 
time series of P^ orus because the deformation of the tori 
near the border makes the preparation of an initial torus 
state difficult [36[. In the present numerical experiment, 
we again take ip torus (q;p c ) as an initial state which is in- 
dependent with f3 (see Appendix |A( with p c = 0.). The 
discrepancy of i(j torus (q; p c ) from an approximate eigen- 
state on the torus produces multiple characteristic time 
scales in P* orus , each of which is associated with an in- 
dividual torus state. Here we estimate 7 from the most 
significant contribution in P* orus [36]. Evidence justify- 
ing such a procedure is shown in Fig. HSIf b): the positions 
of oscillation peaks for the case without noise e = are 
in good agreement with those for the case with small e. 
With increase in e, we can see that sensitive dependence 
of 7 gradually disappears, and 7 becomes larger on av- 
erage. All these are consistent with the results for large 

p. 

Here we show another reasoning that the deformation 
of islands with the changes of s does not induce the sen- 
sitive fluctuations of the amplitude of the tunneling os- 
cillations and 7 in Fig. [15] First, the islands scarcely 
change with noise, as is seen in Fig. [T5] This suggests 
that the influence of the s-dependence of the islands on 
the fluctuations is almost independent of the strength of 
noise. On the other hand, the fluctuation in 7 disap- 
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pears with sufficiently strong noise. Hence we conclude 
that these fluctuations sensitive to s are mainly due to 
the quasienergy resonance. 



VIII. DISCUSSION AND SUMMARY 

In this paper, we have shown that there exists a very 
close relation between dynamical tunneling and dynami- 
cal localization. As mentioned in the introduction, chaos- 
assisted tunneling or chaotic tunneling is usually under- 
stood as a mechanism which enhances tunneling. It is 
certainly true that the transition between quasi-doublets, 
whose support in phase space is symmetric tori, ac- 
quires additional tunneling pathways with the assistance 
of chaotic states. Semiclassical arguments also provide 
a similar picture: exponentially many tunneling paths 
connect the torus region to the chaotic sea. 

On the other hand, what we have seen in the present 
paper is another aspect of the role of chaos in the tunnel- 
ing process. Since the tunneling effect is a typical non- 
local wave phenomenon, it is natural to expect that what 
is happening in chaotic regions must affect the nature of 
the wave function in the integrable region and vice versa. 
It is well known that even without stable invariant com- 
ponents, wave- functions in the chaotic sea are localized 
in nontrivial ways, e.g., scarring due to unstable periodic 
orbits [37]], localization on cantori, and dynamical local- 
ization. One may thus ask the question how and to what 
extent such wave effects appearing in the chaotic region 
affect or are influenced by dynamical tunneling between 
torus and chaos. 

Here we focused particularly on the interplay between 
tunneling and dynamical localization and gave evidence 
that the non-local nature of quantum tunneling becomes 
manifest if one of the quantum effects, that is dynam- 
ical localization, is attenuated. If we add noise to the 
chaotic region, which destroys the interference generat- 
ing dynamical localization, drastic enhancement of the 
tunneling rate is observed. 

The result may be interpreted in two different ways: 
on one hand, we may simply say that chaotic tunneling 
is strongly amplified under the influence of noise. On the 
other hand, we may regard the drastic enhancement of 
the tunneling rate as a result of manifestation of poten- 
tially existing chaotic tunneling. 

The former interpretation is a bit too naive, but most 
arguments concerning chaotic tunneling so far have been 
based on this picture. This implicitly assumes that 
chaotic tunneling is already present even without noise 
and it occurs irrespective of the status of surrounding 
chaotic seas. In this picture, the tunneling process ends 
when the wave function goes through KAM barriers, and 
the behavior shown in Fig. [3] is viewed as manifestation 
of chaotic tunneling. 

In contrast, a semiclassical argument employing com- 
plex classical orbits strongly supports and is consistent 
with the latter interpretation. As mentioned in the in- 



troduction, complex orbits contributing to the tunneling 
amplitude have an amphibious character. They play the 
role of tunneling orbits when they move into the inte- 
grable domain, but behave as almost real orbits after es- 
caping out of the integrable domain. There are exponen- 
tially many complex orbits even in the integrable regime, 
and they flow into the chaotic sea [l4| • Since these orbits 
are almost completely governed by the real dynamics in 
the chaotic sea, the diffusive motion is suppressed by dy- 
namical localization as in the genuine real dynamics. In 
this sense we can say that the tunneling decay is inhib- 
ited due to dynamical localization in the chaotic sea, and 
exponentially many potentially existing tunneling orbits 
appear if noise destroys the coherence creating localiza- 
tion. If this interpretation is true, it should be possible 
to explain the maximal tunneling rate, which was ob- 
served in the strong noise limit (see Fig. [TUJ) and called 
a "classical" tunneling rate, by evaluating the sum over 
complex orbits which flow into the chaotic sea, although 
this evaluation is still not achievable for several reasons 

We further note that, for almost the same reasons, the 
origin of coherent oscillation is not clearly accounted for 
within semiclassical arguments. As presented in Fig. [4j 
the quasi-doublet state was formed between one torus 
state and one chaotic state when scanning system pa- 
rameters. It is reasonable to suppose that some simpler 
form of coherence than destructive interferences generat- 
ing dynamical localization could underlie complex orbits 
connecting torus and chaotic regions, but the mechanism 
invoking the observed strong resonances is not well un- 
derstood. 

The most noteworthy consequence of the interplay be- 
tween dynamical tunneling and dynamical localization is 
that the increase in Planck's constant does not necessar- 
ily lead to the growth of tunneling rate. As Planck's con- 
stant gets larger, the localization length decreases, which 
acts to strengthen the effect of suppression by dynamical 
localization. If the suppression mechanism dominates the 
change in the tunneling rate, a counter-intuitive change 
in the tunneling rate with Planck's constant may happen. 
We verified in Fig. [13] that such a range indeed exists. 

We finally discuss the sharpness of the border between 
torus and chaotic regions. The motivation for construct- 
ing the system so as to have sharply or clearly divided 
phase space was to avoid other complex problems of 
generic phase spaces and make the study of our tunneling 
problem simpler. In particular, the presence of nonlin- 
ear resonances leads to resonance assisted tunneling 
an alternative mechanism enhancing the tunneling rate. 
Conversely, cantori or broken tori play the role of par- 
tial barriers which suppress the transport [H, HO, Slf. It 
seems possible to avoid such problems when we study a 
system with sharply divided phase space. This is cer- 
tainly the case in classical mechanics, but as mentioned 
in section IVHl the sharpness of the boundary, more pre- 
cisely the discontinuity of kinetic or potential functions, 
necessarily invokes diffraction. In section IVIIl we have 
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verified that the claims given for the sharply divided case 
do not nevertheless need to be corrected at least in a 
qualitative level. However, we should recall that the un- 
derlying mechanism is quite different between diffraction 
and tunneling, although both concern classically forbid- 
den processes. We need alternative semiclassical tech- 
niques to obtain quantitative results. In particular, the 
semiclassical evaluation employing complex classical or- 
bits is no longer valid if the system loses analyticity. It 
is also frustrating that the recently proposed formula to 
predict a tunneling rate assumes the extension of tori 
of fictitious integrable system to the surrounding chaotic 
sea, which requires that the border between torus and 
chaotic regions be as sharp as possible [3l[. This alter- 
natively means that the validity of this formula is maxi- 
mally guaranteed in the strong diffraction limit. A quan- 
titative evaluation of tunneling rate is not an easy task 
if we define dynamical tunneling as a process far from 
the diffraction limit. This issue will be dealt with in the 
forthcoming publication [33( ] . 

Throughout this paper, noise is implemented as a 
stochastic process. It is easy to hypothesize that analo- 
gous phenomena will occur when we employ chaotic dy- 
namics to realize noise. Thus we close this paper by 
raising the following question: is the spontaneous recov- 
ery of chaotic tunneling generic in multiple-dimensional 
systems with mixed phase space? 



APPENDIX A: QUANTIZATION OF 
INTEGRABLE TORI 

We show how we approximately construct a torus 
state, which is the initial state of our numerical exper- 
iments (see, Sec. IIII|) . In the torus region, the classical 
dynamics of the system presented in Sec. Ill Al is described 

by 



p n+1 = p n + ksm(2nq n ), 
q n+ i =q n +u. 



(Al) 
(A2) 



We may assume that the momentum of an invariant torus 
is specified by a function of q 



p = P(q). 



(A3) 



Namely, (q,P(q)) is supposed to be on a torus. After 
an iteration of the integrable map, the system arrives at 
(q + u,P(q) + /csin(27rg)), which must agree with (q + 
u, P(q + u)) by the assumption (|A3j) . Hence P(q) must 
satisfy the following equation 



P(q) + k sin(27r<2) = P(q + u) 



(A4) 



To solve it, we expand P(q) as 



P(q) = ^2 ame 



27rimq 



(A5) 



m= — oo 



The functional equation (|A4ft is cast into an algebraic 
equation of a m 



27vimuj 27rimq 



Hence we have 



kcos(2ir(q - |)) 
2 sin(7rw) 



■Pc, 



(A6) 



(A7) 



where an arbitrary constant p c specifies a torus. We re- 
mark that the solution of Eq. (|A4j) exists only when uj is 
not any integer. 

From the Einstein-Brillouin-Keller quantization, we 
obtain the quantized tori 



V tOTUS (<z;Pc) 



= Aexp 
= Aexp 



'2m 



h \47rsin(7rc<;) 



p c q 



(A8) 



where A is a normalization constant. Because 
ip torus (q, p c ) must be single- valued, p c must be quantized. 
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